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We study the thermoelectric coefficients of a multi-level quantum dot (QD) weakly coupled to two 
electron reservoirs in the Coulomb blockade regime. Detailed calculations and analytical expressions 
of the power factor and the figure of merit are presented. We restrict our interest to the limit where 
the energy separation between successive energy levels is much larger than the thermal energy (i.e., 
the quantum limit) and we report a giant enhancement of the figure of merit due to the violation of 
the Wiedemann- Franz law when phonons are frozen. We point out the similarity of the electronic 
and the phonon contribution to the thermal conductance for zero dimensional electrons and phonons. 
Both contributions show an activated behavior. Our findings suggest that the control of the electron 
and phonon confinement effects can lead to nanostructures with improved thermoelectric properties. 
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I. INTRODUCTION 

The thermoelectric properties of nanomaterials have 
attracted a great deal of interest in the last few years 
due to the potential applications in power generation and 
refrigeration There is increasing evidence that the 
thermoelectric performance of nanostructured materials 
can be significantly improved. The efficiency of a device 
to convert heat into electricity and vise- versa is measured 
by the dimensionless figure of merit ZT given by 



ZT 



S 2 GT 



(1) 



where S is the thermopower, G is the electrical conduc- 
tance, T is the absolute temperature and n is the thermal 
conductance. There are two contributions to n: the elec- 
tronic, K e , and the phonon, n p h- 

Good thermoelectric materials are considered to be 
those with ZT > 3 at room temperature 0]. The best 
commercial thermoelectric material is bulk Bi2Te3 and 
its alloys with Sb and Se with ZT close to 1. Although 
there is no thermodynamic requirement that imposes an 
upper limit on ZT, until 2001 no progress had been made 
in increasing ZT above 1. The difficulty is due to the in- 
terrelations between the coefficients S, G and k that set 
limitations in varying them independently. A character- 
istic example is the Wiedemann-Franz law according to 
which the ratio n e /GT remains constant. 

Since 2000 there has been a number of experiments in 
nanostructures where the barrier of ZT = 1 has been 
overcome at high temperatures. Venkatasubramanian et 
al. measured ZT « 2.4 in p-type Bi 2 Te3/Sb 2 Te3 su- 
perlattices at T = 300 K. The reason for the enhanced 
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figure of merit is the phonon blocking due to acoustic 
mismatch between Bi2Te3 and Sb2Te3 that causes a sig- 
nificant reduction to n p h- Harman et al. [4| reported 
values of ZT between 1.3 to 1.6 in PbSeTe/PbTe quan- 
tum dot superlattices. More recently Kanatzidis and 
co-workers found that bulk AgPbisSbTe2o with inter- 
nal nanostructures has ZT « 2 at T = 800 K @. In 
nanocrystalline BiSbTe bulk alloys ZT reached the value 
1.4 at T = 373 K Q . The enhanced ZT was attributed to 
the significant reduction of K p h due to the strong phonon 
scattering by the interfaces between the nanostructures. 
A 100-fold improvement of ZT compared to the bulk 
value has been reported recently in Si nanowires 0, H| ■ 
In [7| the increased ZT is due to particularly small values 
of Kph that reach the amorphous limit, while in [8] the 
phonon-drag effect 0, [l(| was identified as the reason for 
the increased thermopower at T = 200 K. 

Early theoretical work pointed out the possibility of 
enhancing ZT in structures of reduced dimensionality 
due to the increased density of states near the Fermi 
level [ll[. Mahan and Sofo jl2| predicted a maximiza- 
tion of the thermoelectric efficiency in materials with a 
delta-like density of states. This suggests that quantum 
dots and molecular junctions are promising candidates 
for good thermoelectric materials. Humphrey and Linke 
fl3j provided a thermodynamic explanation why the op- 
timum density of states for enhanced ZT is a delta func- 
tion. More recent theoretical studies [l4T - [T8l | point out 
the possibility of significantly improving the thermoelec- 
tric efficiency in small quantum dots and molecules. 

The present work presents a comprehensive analy- 
sis of the thermoelectric properties of a weakly coupled 
multi-level quantum dot (QD) in the sequential tunneling 
transport regime and suggests the possibility of a huge in- 
crease of ZT in the Coulomb blockade regime due to the 
violation of the Wiedemann-Franz law. We focus on the 
quantum limit where the energy separation between sue- 
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cessive electronic levels is large compared to the thermal 
energy. Preliminary results have been presented in a con- 
ference form [IE] . Deviations from the Wiedemann- Franz 
law in zero dimensional structures were predicted also in 
[HlHoI without, however, relating these to the possibility 
of an enhancement of ZT. Recently, an interesting pre- 
diction for a significant increase of ZT in a single-level 
molecule was made by Murphy et al. EH . In Ref. [l6| it 
was found that the energy parameter that controls ZT is 
the electron-electron repulsion. Here we show that for a 
multi-level system in the Coulomb blockade regime ZT 
is controlled by the level spacing and, consequently, the 
spatial confinement. 

The Coulomb blockade oscillations of the conductance 
and thcrmopower in the sequential tunneling regime have 



been studied in the seminal papers by Beenakker 21 1 
an d by Beenakker and Staring[22]. Recently, the theory 
in [2l[ |22T ] , that is based on a master equation approach 
within the constant interaction model, has been extended 
to calculate the electronic contribution to thermal con- 
ductance [3 EE Ell- Here we present, for the first time, 
a detailed derivation for an analytical expression for re e 
which serves as a solid framework for the understanding 
of the behavior of this coefficient in the quantum limit. 
Within this framework we discuss the validity of previous 
results [23( 1 . The cases of a QD with equidistant and 
non-equidistant energy spectrum are examined. We also 
present detailed calculations and analytical derivations 
for the power factor S 2 G and the figure of merit. We 
predict an optimization of the thermoelectric efficiency 
at specific values of the Fermi level. An important out- 
come of this work is the similarity in the temperature 
dependence of the electronic and phonon contribution to 
k in structures with zero dimensional (OD) electrons and 
phonons. Our findings suggest the possibility of a dra- 
matic increase in ZT in both electron- and phonon- block- 
ing devices. 



of the heat current Q through the dot when small temper- 
ature and voltage differences AT and AV, respectively, 
are applied between the two reservoirs. The linear re- 
sponse model proposed in the pioneer work by Beenakker 
21] and Beenakker and Staring [22| for the calculation 
of the charge current J has been generalized to include 
Q 0, EE El- Then Q is given by 



L-nR 



Q = - 



si 



eAV 
k^T 



AT \ 



(2) 



In the above equation the energy dependence of the tun- 
nel rates has been ignored, e is the magnitude of the 
electron charge. Moreover, for the quantities s m with 
to = 0, 1, 2 we use the following notation 



00 00 



s m = J2J2 Pe q (N)P eq (e n \N)[l-f(E-E F )}(E-E F ) m 

n=l N=l 

where, the sum indexes n and N refer to the energy level 
and the number of electrons in the dot, respectively, e n is 
the energy of the n state and E — e n + U(N) — U(N—1) 
with U(N) being the electrostatic energy for a dot with N 
electrons. Moreover, P eq (N) is the probability that the 
dot contains N electrons in equilibrium and P e q(e n \N) is 
the probability that in equilibrium the n level is occupied 
given that the dot has ./V electrons. The expressions for 
P eq (N) and \N) are given in [22j. 

Utilizing the standard transport equation (Tol j 



= MAV + KAT 



(4) 



the transport coefficients M and K can be readily ob- 
tained pja], namely, 



M 



kuT 



T tot Sl 



(5) 



and 



II. THEORY 

We consider a quantum dot that consists of non- 
degenerate discrete energy levels e„ (n = 1,2,3,...) 
weakly coupled to two electron reservoirs. In equilib- 
rium both reservoirs have the same temperature T and 
the same chemical potential Ep. The energy distribution 
in the reservoirs is the Fermi-Dirac function f(E — Ep) = 
l/{exp[/3(£ - E F )\ + 1}, where /3 = (k B T)-\ We ex- 
amine the Coulomb blockade regime where the charging 
energy Eq = e 2 /2C (C is the capacitance of the dot to 
the surroundings) is larger than the energy spacing Ae 
between successive energy levels. The tunnel rate from 
level n to the left, L, (right, R,) reservoir is denoted as 
r« (r^)- We assume that ksT and Ae are much larger 
than h(Tn + (where h is Planck's constant) and we 
neglect effects due the virtual tunneling (cotunneling) of 
electrons through the dot [2l|, [22[ . 

The results presented here are based on the calculation 



K 



k B T 2 



T tot s 2 . 



(6) 



In the above equations T tot = T L T R /{T L + T R ). 

Onsager's symmetry relation relates the thermoelectric 
coefficient M to the thermopower S by M=SGT 
Then S is written as 



s 



M 



1 sx 



GT eT s ' 



where G is given by [2l| 



G — Gqsq, 



(7) 



(8) 



with G = (e 2 /k B T)T tot . The expression we derive 
for S is identical to the expression (3.13) obtained by 
Beenakker and Staring [22[ , where the authors calculated 
the electron current through the dot instead of the heat 
current. The approach we use here for S is the so-called 
II-approach fioj ]. 
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Finally, the electronic contribution to the thermal con- 
ductance can be calculated by using the standard rela- 
tion o 



= -K- S 2 GT. 



(9) 



In Fig. la we show the calculations of the thermoelec- 
tric coefficient M as a function of the Fermi level for 
a three- level QD when (3Ae = 10 and Ae = 0.5-E C - 
M exhibits Coulomb blockade oscillations with period- 
icity AEf — Ae + (e 2 /C) as does the electrical con- 
ductance [2l|. The effect of temperature on the shape 
of M is shown in Fig. lb. The Coulomb blockade os- 
cillations of the coefficient K are shown in Fig. 2a for 
the same parameters. As we can see —K shows a dou- 
ble peak structure which is reproducible in energy in- 
tervals Ae + (e 2 /C). The theoretical values for K have 
been divided by L TG max where L = (k B /e) 2 n 2 /3 and 
Qmax _ (2 Q /4 i s t ne maximum value of the electrical con- 
ductance. The effect of temperature on the peaks of — K 
is shown in Fig. 2b and it will be explained by a simple 
analytical expression derived later. 

In Fig. 2 we have also shown the calculation of S 2 GT. 
The curves for —K and S 2 GT are indistinguishable so 
that At e is near zero. It is found that the magnitude of 
the maximum values of n e collapses rapidly to zero as 
(/3Ae) 2 exp(— /3Ae) violating the Wiedemann-Franz law. 
Usually, in degenerate 2D and ID systems S 2 GT is a 
small quantity compared to —K and the thermal con- 
ductance follows closely —K. The standard Sommerfeld 
expansion leads to the Wiedemann-Franz law according 
to which K e = —K = LqGT . The above imposes a sig- 
nificant limitation on ZT when S 2 < Lq even when the 
phonon contribution to the thermal conductance is ig- 
nored. 

The similarity between — K and S 2 GT can be ex- 
plained in simple terms. Namely, the dominant con- 
tribution to the summations over n and N in Eq. ^ 
occurs for n — N — N m i n where N m i n is the integer 
which minimizes the absolute value of the energy differ- 
ence A(N) = e N + U(N)- U{N -l)-E F .We note that 
A(N m i n ) oscillates in a sawtooth manner with Ep and 
passes through zero each time an electron enters the dot. 
In what follows we keep Beenakker's notation [2l[ and 
we write A(N m i n ) = A m i n . The oscillations of A m i n are 
shown in Fig. 3. 

In the low-T limit the probabilities P ea (N m j n ) and 
Peq{f-N min \N m in) can be approximated by |2l| 

PeqiNmvn) « /(A mi „) (10) 



and 



Peq{e-N mm \N min ) W 1 



(11) 



Inspection of Eq. Q shows that the quantities s m can 
be written in the approximate form 
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FIG. 1: The Coulomb blockade oscillations of M and the 
effect of temperature on the shape of M for a three-level 
dot. The solid line in (a) is for k B T = 0.05E C . In (b) 
The solid, dashed-dotted and dotted lines correspond to 
k B T = 0.05,0.035 and 0.02E C , respectively. The level spac- 
ing is Ae = 0.5Ec- 



where the superscript N m i n denotes the contribution of 
the term n = N m i n in the sum ([3]). Now, by substituting 
Eq. (TT2"]) into Eqs.© and © one readily gets 



S 



M 
GT 



1 

"eT" 



-,A T 



(13) 



and 



S 2 GT = -K 



k B T 2 



(14) 

We recall that f(A min ) = l/[exp(/3A mm ) + 1]. 

By maximizing Eq. (114[) with respect to A m „ we find 
that —K and S 2 GT exhibit two symmetric peaks at 
(3A m i n = ±2.4 shown in Fig. 2. The magnitude of the 
peaks is 



(S 2 GT) max = 0.44/fc B r' 



(15) 



f \Afyiin) [1 f (A m i n ^~\ (A m ^ n ) , (12) 



We note that the position of the two strong maxima are 
in agreement with the predicted maximization of S 2 GT 
and ZT in materials with delta-like density of states [12j . 
Moreover, according to Eq. (fl~3"|) . the absolute magnitude 
of S that corresponds to the optimization of S 2 GT is 
2Aks/ e = 207 /iV/K which is in absolute agreement with 
the value predicted in (l2l |. 

At this point we should remark that Eq. (f]"2"|) describes 
accurately the magnitude and the shape of the physical 
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FIG. 2: The coefficient K as a function of Ef for a three-level 
dot. The calculated values of K are divided by L TG max 
where Q max = Go/4 is the maximum value of the conduc- 
tance. In (a) the Coulomb blockade oscillations of K are 
shown for ksT = 0.05-Ec- The effect of temperature on the 
shape of K is shown in (b). The solid, dashed-dotted and 
dotted lines correspond to ksT = 0.05, 0.035 and 0.02Ec, re- 
spectively. The level spacing is Ae = 0.5-Ec. In the same 
figure we plot S 2 GT. The curves for -K and S 2 GT are in- 
distinguishable. 
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FIG. 3: The Coulomb blockade oscillations of A m i n . Ae = 
0.5E C . 



quantities G, M, K and S 2 GT. Moreover, the expres- 
sion (ITUl) for S is in very good agreement with the nu- 
merical result around A m £„ = within the energy range 
I ^min | ^ AE. The contribution of the terms n ^ N m i n 
in Eq. © is responsible for the fine structure on the oscil- 
lations of S that was predicted by Beenakker and Staring 



22]. These terms do not affect S 2 GT and are not dis- 
cussed here for simplicity. 

The strong competition between —K and S 2 GT in 
Eq.© gives n e = if only the leading contribution to 
the sum ([3]) for n = N m i n and N = N m i n is considered. 
The terms that correspond to n = N m i n ± 1 disturb the 
above cancellation and, as we show below, an analytical 
expression for n e can be obtained. Inspection of Eqs.© 
and (|10p shows that these terms introduce tiny correc- 
tions to the quantities s m of the form 



,±1 



,±1 



= P eq {e N . 
x /(A mm )[l 



\N mm )(A mm ± Ae) r 
f(A mm ±Ae)} 



(16) 



where the superscripts AT min ± 1 denote the contribu- 
tion from the terms n — N m i n ± 1. The probabilities 
Peq{tN mm ±i\N mm ) are given by [22| 



Peq{^N min +i\N min ) w exp(-/3Ae) 



and 



Pe q {e Nmm -i\N n 



1. 



(17) 



(18) 



In the energy interval around A m i n = 0, which is 
the regime where G exhibits the well known peaks [2l| . 
s m mi ' l±1 are smaller by a factor of exp(—/3Ae) compared 
to s„""™ given by Eq. (O (we note that /3Ae > 1). 
This explains why these corrections are unimportant for 
G [2l|. To obtain n e we use Eq.© after we substitute 
K, S and G from Eqs.@, © and ©, respectively. Then 
we take 



r' 



k R T 2 



■S'2 



so 



(19) 



We, now, approximate s m by the sum s^ mi -r o m 
Sm min+1 an d we handle s^ m " 1=tl as small quantities in 
order to linearize expression (|19p by keeping terms up to 
first order in exp(— /3Ae) and obtain 



where, 



and 

K N min ±l 



k B T 2 



b 2 



N min 
S 



(20) 



(21) 



M mi „±l 



k B T 21 



(4 



N min N mirl ±l 
s 



2s 



N„ 



JV m( „±l 



)]■ 

(22) 



The substitution of the quantities s^" 
from Eqs. (fT2")) and (p~^]) gives 



and s 



s N mln ±i = k B T tot P eq (e Nmm±1 \N mm )(/3Ae) 2 
x /(A mm )[l-/(A mm ±Ae)] 



(23) 



5 




In deriving the above equation we have taken into ac- 
count that exp(/3Ae) 3> 1 in the quantum limit. Eq. (|27|) 
denotes a strong sensitivity of n e to the details of the 
energy spectrum. 

Inspection of Eqs.® and (|27|) shows that the 
Wiedemann-Franz law is not valid in the quantum 
limit. We find that the peak values of the Coulomb 
blockade oscillations of smaller than that the 

Wiedemann-Franz law predicts by a factor of the order 
(/3Ae) 2 exp(— f3Ae). Namely, the ratio of the peak values 
of K e and G for A m ;„ = is given by 
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FIG. 4: The two symmetric contributions to K e (solid line) 
for a dot with equidistant energy spectrum around Amin = 
0. The dashed-dotted and the dotted lines correspond to 
K N min -l an( j K N min +i^ reS p ec tively. Ae = 0.5Ec and 
/3Ae = 15. ko = L TG max (pAe) 2 exp(-/3Ae) where G max = 
Go /4. The shown structure is reproducible in energy intervals 
Ae + (e 2 /C) except for the first and the last level where « e is 



_1 and ftp 



respectively. 



It is apparent from the above analysis that n e consists of 
two contributions 



Mmin-1 



(24) 



These contributions are shown as dashed-dotted and dot- 
ted lines in Fig. 4. We see that they are symmetric around 
A m i„ = and peak at A min = — Ae/2 and Ae/2. We 
note that K^ min+1 and re^" 1 "* -1 are periodic as a function 
of Ep with periodicity Ae + (e 2 /C). 

We should mention that there is a symmetry in the 
correction terms G Nmin±1 for the electrical conductance, 



arising from the terms n = N„ 
Namely 



± 1, and 



,±1 



— ] T((3Ae) 2 G 
e 



N„ 



,±1 



,±1 



(25) 



Such a symmetry was absent in the papers [la , |23| | . This 
poses a question about the existence of a solid theoretical 
background for the expressions for n e that appeared in 
these papers and about the validity of the relevant results 
especially for the case of a dot with degenerate energy 
spectrum. 

Substituting Eq . (|23|) in (|24l) . after some trivial algebra 
we take tH, fijl 



k B T tot (pAe) 2 
(H 



3 /3A e ( e 2/3A mI „ + X ) + e 0A min ^20Ae + 



.(26) 



Eq. (|26p is in excellent agreement with the numerical 
result. The magnitude of the peak values of n e around 
Amin — is given by the simple expression 



k B r tot (/3Ae) 2 exp(-/3Ae). 



(27) 



G" 



L T(^Ae) 2 exp(-/3Ae). 



(28) 



As we show below the violation of the Wiedemann-Franz 
law in the system under consideration leads to an ex- 
ponential increase of ZT with /3Ae. Also, recently it 
has been reported that deviations from the Wiedemann- 
Franz law lead to enhancement of ZT in nanocontacts 
made of two-capped single wall nanotubes [2 lj , quantum 
dots attached to ferromagnetic leads [HI], and strongly 
correlated quantum dots[25j. 

The maximization of S 2 GT , discussed earlier, has a di- 
rect impact on the dimensionless figure of merit ZT when 
the phonon contribution to the thermal conductance is 
ignored 



ZT = 



S 2 GT 



(29) 



According to the above expression and Eqs.([T5j) and (|27|) 
ZT shows two strong maxima at f3A m i n — ±2.4 of mag- 
nitude 



(ZT) r 



{S 2 GT) r 



= 0.44exp(/3Ae)/(/3Ae) 5 



(30) 

The structure of ZT is shown in Fig. 5c. The exponential 
dependence of these maxima can be readily explained by 
the fact that, in the energy interval where S 2 GT exhibits 
the peak structure shown in Fig. 5a, K e remains constant 
and its magnitude is given by Eq. (|27| (see Fig. 5b). 

Eq. (30) shows that the thermoelectric efficiency can 
be dramatically enhanced by controlling the size of the 
dot and consequently Ae. For small dots (or molecules) 
the requirement that the system is in the quantum limit 
is fulfilled even at room temperature. Practically speak- 
ing, the quantum limit is approached when f3Ae > 5. 
High values of ZT in the range of 2.6 to 97 at T = 300 K 
(fcgT = 0.026 eV) can be obtained when Ae varies be- 
tween 5 to 10 k B T (namely, between 0.13 to 0.26 eV). We 
note that the energy separation between the first two en- 
ergy levels in a spherical dot for a material with electron 
effective mass of the order of 0.1 m e (m e is the electron 
mass) is greater than 0.15 eV when the diameter is less 
than 5 nm. Finally, we add that in our model we have 
assumed that the charging energy Eq is larger than Ae. 
Large values of Ec of the order of the tenth of eV have 
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FIG. 5: Calculated values of S 2 GT, K e and ZT as a func- 
tion of Ef for a three-level dot. The level spacing is 
Ae = 0.5i?c. The solid and dashed lines correspond, re- 
spectively, to /3Ae=ll and 9. The results are plotted around 
Ef = 4Ec that corresponds to A m in = when an elec- 
tron enters the dot at the level n = 2. In (b) ko = 
L TG maw (PAef exp(-^Ae) where G maa: = G /4 and in (c) 
(ZT) max = 0.44exp(/3Ae)/(/3Ae) 2 . 



been reported recently in single-electron transistors with 
organic molecules acting as Coulomb islands (23 l27|. 

In a recent paper [23[ the possibility of tuning ZT of 
a dot by varying was recognized. In [23| no detailed 
calculations were performed and the possibility of tun- 
ing ZT was incorrectly related solely to the dependence 
of K e on the energy spectrum and the variation of n e 
with Ep. However, the variation of S 2 GT with Ep is 
the most remarkable behavior predicted in these systems. 
Although, the phonon contribution to the thermal con- 
ductance could washout the variation of the total thermal 
conductance as a function of Ep, ZT will still vary peri- 
odically with Ep following the structure of S 2 GT shown 
in Figs. 2 and 5a. Moreover, the maximization of S 2 GT 
could be utilized for the prediction of an optimal value 
for K p h that maximizes ZT [16j |, 

So far our analysis has referred to a dot with equidis- 
tant energy levels. The generalization to a non- 
equidistant energy spectrum is straightforward. How- 
ever, we should remark that this generalization will not 
affect our findings for S 2 GT (apart from the absence of a 
regular periodicity in the oscillations of S 2 GT) but it will 
only affect the shape of ZT in the case where the phonon 
contribution to thermal conductance is small compared 



FIG. 6: Calculated values of S 2 GT, n e and ZT as a func- 
tion of Ef for a dot with non-equidistant energy spectrum. 
The results are plotted around -Ef = 4_Ec that corresponds 
to A m i n = when an electron enters the dot at the level 
n = 2. The energy separation between the two adjacent 



levels is Ae B 



0.5E C and Ae ne = 0.55E C . /3Ae p 
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and /3Ae ne = 16.5. The dashed-dotted and the dotted 
lines in (b) show the contributions n^ min ~ 1 and n^ min+1 , 
respectively. k = L TG max (/3Ae ne ) 2 exp(-/3Ae ne ) and 
{ZT) max = 0.44exp(^Ae ne )/(^Ae„ e ) 2 . 



to n e . This is due to the profound effect of the energy 
spectrum on the shape of n e . In Fig. 6 we present calcu- 
lations for S 2 GT, K e and ZT for a three-level dot when 
one electron enters the dot at the level n — 2. The en- 
ergy separation between the next and the previous energy 
level is Ae„ e — 0.55Ec and Ae pr — 0.5Ec- To obtain n e 
in the case of dot with non-equidistant energy spectrum 
one needs to replace Ae by Ae ne and Ae pr in the expres- 
sions (|2"5)) for K^ min+1 and K^ mi " _1 , respectively. In this 
case the two contributions to n e become asymmetric as 
shown by the dashed-dotted and dotted lines in Fig. 6b. 
This results to the two asymmetric peaks in ZT shown 
in Fig. 6c. 

We should remark at this point that the phonon con- 
tribution to thermal conductance could significantly sup- 
press ZT . However, the phonon confinement in free- 
standing QDs or QDs embedded in a matrix material, 
where there is a significant mismatch in the elastic prop- 
erties of the QD and the surrounding material, could dra- 
matically decrease the phonon contribution K p h- The 3D 
confinement of phonons in nano-objects results in sim- 
ilar restrictions in the phase space of the phonon wave 
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vector as in the case of electrons. The discrete phonon 
frequency spectrum influences the phonon heat transfer 
through a nanostructure. In what follows we will attempt 
to reveal the similarities between the fermionic and the 
bosonic contributions to the thermal conductance. As we 
will show both contributions show an activated behavior. 

We start from the Landauer expression used by Rego 
and Kirczenow 28] for the heat transferred by ID bal- 
listic phonons through a quantum wire attached to two 
thermal reservoirs in the presence of a small tempera- 
ture difference AT between them. This is modified here 
to account for the discrete phonon spectrum of the QD. 
We assume that the temperature in the right reservoir 
is raised by AT compared to the temperature in the left 
reservoir. Then the phonon contribution to the heat flow 
from the right to the left reservoir is written in the form 

m 



Qph 



E 



dojhujTi(uj)[nR(uj) — riL{u>)] 



(31) 



where T(w) denotes the transmission coefficient for a 
phonon of frequency w through the ith-mode into the 
nanostructure and tlr{uj) = {^p[ kB ^+AT} ^ ~ l} -1 anc ^ 
jil(w) = [exp(^y) — are the phonon distribution 
functions in the right and the left reservoirs, respectively. 
The summation is over all the vibrational modes. For the 
case of resonance energy transfer we write the transmis- 
sion coefficient T(cj) in the form 30] 



TH 



r 



pi y r 



(32) 



ph ph 

where T p h denotes the phonon tunnel rate. Then we take 



(33) 



where r$ = I*r« /(I* + r« ). The expression © 



ph phi \ ph 



ph! 



is in agreement with the expression derived recently by 
Segal and Nitzan [3(| for the heat current through a har- 
monic molecule weakly coupled to two thermal reservoirs 
by using a master equation approach. 

In the linear regime, where AT is small, it can be easily 
shown that 



dn L (ui)AT 
d{fiu;i) T ' 



where. 



-Pi 



d{hjji) [exp(/3fiwi) — l] 2 

At temperatures where flftidi >> 1 we can write 
dn L (wi) 



d{f\u)i) 



-0exp(-/3tkJi) 



(34) 



(35) 



(36) 



Then, the phonon contribution to the thermal conduc- 
tance K p h = Qph I AT is written in the form 



x ph 



= k B J2 r^CSfiWi) 2 expi-ptkJi). 



(37) 



At low T only the first low energy phonon mode Huiq 
contributes significantly to K p h- Then we can write 



K P h = kBTffififuuo) 2 exp(-phiuo). 



(38) 



The similarity of the above expression with Eq. (|2"7} for 
K e is quite remarkable. Interestingly, as in the case of 
the ID thermal conductance [28|, the phononic and the 
electronic contributions to the thermal conductance in 
structures with a 3D electron and phonon confinement 
appear to show a similar T-behavior. 

According to our analysis the maximum value for ZT 

is 



(ZT) r 



0.44fc B r*°' 



(39) 



where n e and n p h are given by Eqs. ([2"Tj) and (|3"5j) . Due 
to the activated terms in these equations ZT can take 
values much larger than 1. Our finding suggests that the 
control of the electron and phonon confinement effects in 
nanostructures can improve significantly their thermo- 
electric properties. Namely, for tkuo and Ae of the order 
of 150 to 200 meV and r tot « Y p f we find that ZT can 
reach the values 2.2 to 8.4 at T = 300 K. We also note 
that (ZTy nax is an increasing function of r to< 



t lytot 
I L ph- 



III. CONCLUSIONS 

In summary, we have presented a comprehensive study 
of the thermoelectric coefficients of a multi-level QD 
weakly coupled to two electron reservoirs in the Coulomb 
blockade regime for sequential tunneling. Our analysis 
focuses on the quantum limit (2l| where the level spac- 
ing between successive electronic levels is much larger 
than the thermal energy (/3Ae >> 1). Analytical ex- 
pressions for the power factor and the figure of merit 
are derived. We found an exponential increase of ZT 
with /3Ae when phonons are ignored. This is due to 
the activated behavior of K e . We also show that n p h 
for a dot with discrete phonon spectrum shows a simi- 
lar activated behavior. Interestingly both fermionic and 
bosonic contributions to the thermal conductance show a 
very similar T-dependence. Although, it was previously 
believed that electron-transmitting and phonon-blocking 
structures have good thermoelectric properties our 
results now open a route for designing electron-blocking 
and phonon-blocking nanostructures with improved ther- 
moelectric performance. 



8 



[1] M. S. Dresselhaus, G. Chen, M. Y. Tang, R. G. Yang, H. 

Lee, D. Z. Wang, Z. F. Ren, J.-R Fleurial, and P. Gogna, [14 
Adv. Mater. 19, 1043 (2007). 

[2] A. Majumdar, Science 303, 777 (2004). [15 

[3] R. Venkatasubramanian, E. Siivola, T. Colpitts, and B. 
O'Quinn, Nature 413, 597 (2001). 

[4] T. C. Harman, P. J. Taylor, M. P. Walsh, and B. E. 

LaForge, Science 297, 2229 (2002). [16 

[5] K. F. Hsu, S. Loo, F. Guo, W. Chen, J. S. Dyck, C. Uher, 

T. Hogan, E. K. Polychroniadis, and M. G. Kanatzidis, [17 
Science 303, 818 (2004). 

[6] B. Poudel, Q. Hao, Y. Ma, Y. Lan, A. Minnich, B. Yu, [18 
X. Yan, D. Wang, A. Muto, D. Vashaee, X. Chen, J. Liu, 
M. S. Dresselhaus, G. Chen, and Z. Ren, Science 320, [19 
634 (2008). [20 

[7] A. I. Hochbaum, R. Chen, R. D. Delgado, W. Liang, E. 

C. Garnett, M. Najarian, A. Majumdar, and P. Yang, [21 
Nature 451, 163 (2008). [22 

[8] A. I. Boukai, Y. Bunimovich, J. Tahir-Kheli, J. -K. Yu, 

W. A. Goddard III, and J. R. Heath, Nature 451, 168 [23 
(2008). [24 

[9] R. Fletcher, E. Zaremba, and U. Zeitler in Electron- 

Phonon Interactions in Low Dimensional Structures, [25 
edited by L. Challis (Oxford University Press, Oxford, 
2003), p. 149. [26 
[10] M. Tsaousidou in The Oxford Handbook of Nanoscience 

and Technology, edited by A.V. Narlikar and YY. Fu, [27 
Vol.11, Ch.13, (Oxford University Press, Oxford, 2010), 
P-477. [28 
[11] L. D. Hicks and M. S. Dresselhaus, Phys. Rev. B 47, 

12727 (1993); L. D. Hicks and M. S. Dresselhaus, Phys. [29 
Rev. B 47, 16631 (1993). 
[12] G. D. Mahan and J. O. Sofo, Proc. Natl. Acad. Sci. USA [30 

93, 7436 (1996). 
[13] T. E. Humphrey and H. Linke, Phys. Rev. Lett. 94, 



096601 (2005). 

M. Tsaousidou and G. P. Triberis, 
|arXiv:cond-mat /0605286"] (May 2006). 
M. Tsaousidou and G. P. Triberis, in Physics of Semicon- 
ductors: 28th International Conference on the Physics of 
Semiconductors-ICPS 2006, AIP Conf. Proc. 893 (AIP, 
New York 2007) p. 801-802. 

P. Murphy, S. Mukerjee, and J. Moore, Phys. Rev. B 78, 
161406(R) (2008). 

C. M. Finch, V. M. Garcfa-Suarez, and C. J. Lambert, 
Phys. Rev. B 79, 033405 (2009). 

R. Swirkowicz, M. Wierzbicki, and J. Barnas, Phys. Rev. 
B 80, 195409 (2009). 

X. Zianni, Phys. Rev. B 75, 045344 (2007). 

B. Kubala, J. Konig, and J. Pekola, Phys. Rev. Lett. 
100, 066801 (2008). 

C. W. J. Beenakker, Phys. Rev. B 44, 1646 (1991). 

C. W. J. Beenakker and A. A. M. Staring, Phys. Rev. B 
46, 9667 (1992). 

X. Zianni, Phys. Rev. B 78, 165327 (2008). 

K. Esfarjani, M. Zebarjadi, and Y. Kawazoe, Phys. Rev. 

B 73, 085406 (2006). 

T. A. Costi and V. Zlatic, Phys. Rev. B 81, 235127 
(2010). 

Y. Wakayama, T. Kubota, H. Suzuki, T. Kamikado, and 

S. Mashiko, Journal of Appl. Phys. 94, 4711 (2003). 

Y. Wakayama, T. Kubota, H. Suzuki, T. Kamikado, and 

S. Mashiko, Nanotechnology 15, 1446 (2004). 

L. G. C. Rego and G. Kirczenow, Phys. Rev. Lett. 81, 

232 (1998). 

A. Ozpineci and S. Ciraci, Phys. Rev. B 63, 125415 
(2001). 

D. Segal and A. Nitzan, Phys. Rev. Lett. 94, 034301 
(2005). 



